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Abstract. Dielectric relaxation in disordered dielectric mixtures are presented 
by emphasizing the interfacial polarization. The obtained results coincide with 
and cause confusion with those of the low frequency dispersion behavior. The 
considered systems are composed of two phases on two-dimensional square and 
triangular topological networks. We use the finite element method to calculate 
the effective dielectric permittivities of randomly generated structures. The 
dielectric relaxation phenomena together with the dielectric permittivity values at 
constant frequencies are investigated, and significant differences of the square and 
triangular topologies are observed. The frequency dependent properties of some of 
the generated structures are examined. We conclude that the topological disorder 
may lead to the normal or anomalous low frequency dispersion if the electrical 
properties of the phases are chosen properly, such that for "slightly" reciprocal 
mixture--when o~i 2> <T2, and e± < e2-normal, and while for "extreme" reciprocal 
mixture-when a\ 3> <T2, and t\ ^-anomalous low frequency dispersions are 
obtained. Finally, comparison with experimental data indicates that one can 
obtain valuable information from simulations when the material properties of the 
constituents are not available and of importance. 



1. Introduction 

Numerical simulations and analytical formulas on dielectric behavior of materials with 
two different phases (binary mixtures) on regular two dimensional structures have 
been shown to exhibit dielectric behavior of Debye type when the concentration, q, 
of inclusion phase is low (Tuncer et al. 2001a, Tuncer et al. 2002b). The polarization 
relaxation appears due to the differences between the electrical properties of the 
constituents, t\l<J\ 7^ £2/02- This polarization is called Maxwell- Wagner-Sillars 
(mws) (Maxwell 1891, Wagner 1914, Sillars 1937) polarization (also known as the 
interfacial polarization) . At higher concentrations of inclusions the dielectric behavior 
deviates slightly from Debye type (Tuncer et al. 2001a) when the mixture is a match- 
composite-the electrical properties of phases, permittivity e and direct-current (dc) 
conductivity c, are chosen in such a way that s\ < £2 and a\ < 02 where the subscripts 
(1) and (2) indicate the matrix and inclusion phases. The deviation from Debye type 
increases for reciprocal-composites-when e\ < £2 and a\ > 02 (Tuncer et al. 2002a). 
Previously, mws polarization in disordered structures have been reported (Tuncer et 
al. 2002a), and non-Debye dielectric relaxations have been observed. 
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In this paper, we stress on the dielectric properties of disordered binary 
mixtures, and extend the investigations on mixtures to micro-structural differences and 
importance of ratio of electrical properties of constituents. The electrical properties of 
the structures were calculated by using the finite element method (fem) (Tuncer 2003). 
The dielectric responses of several selected disordered systems are presented. The 
dielectric permittivity values calculated were as high as in the porous solids (Nettelblad 
1996, Papathanassiou and Grammatikakis 2000, Papathanassiou 2001, Capaccioli et 
al. 2001, Shen et al. 1990). Therefore, the results of two dimensional calculations 
are compared not only to dielectric responses of binary polymeric composites (Tuncer 
and Gubahski 2000, Tuncer and Gubahski 2001b) and of porous liquid-solid systems 
as well. 

2. Low frequency dispersion 

Dielectric properties of materials can be characterized by the complex dielectric 
permittivity, 

e{u) = e'{uo) - ie"{uo) 

or in susceptibility notation, 

x(uj) = e{uj) - e 

which have the angular frequency to dependencies. Here, the instantaneous dielectric 
permittivity e consists of the relaxations which are faster than the relaxations observed 
in the considered frequency range. External variables, such as, temperature, pressure, 
humidity etc., also influence the dielectric properties of materials. When a relaxation 
process is present, the dielectric response of materials usually shows a peak in the 
frequency dependence of the imaginary part of the loss component, x"( w )j and a 
dielectric increment in the real component, x'( w ) (Jonscher 1983). The response is 
often of non-Debye type (Ramakarishnan and Lakshmi 1987), and can in general be 
expressed by Havriliak-Negami empirical formula (Havriliak and Negami 1966), 

e{u;e, X (P),T,a,ff) = e+ M f°\ aW (1) 

[1 + (lUJT) a \ P 

where r and x(0) represent the average relaxation time (inverse relaxation rate) and 
the dielectric strength of the polarization, respectively, a and j3 are parameters that 
display the degree of deviation from the Debye relaxation, a < 1 A a[3 < 1. When 
a = (3 = 1, Eq. JTJ yields Debye's relaxation equation (Debye 1945). The other 
empirical expressions that can be obtained from Eq. (JIJ are when (3=1, Cole- 
Cole (Cole and Cole 1941) and when a = 1 Davidson-Cole (Davidson and Cole 1951), 
respectively. 

At frequencies much lower and much higher than the inverse relaxation time, r _1 , 
the susceptibility can be written as 

X'(w) oc x"H oc w~ a P w»r- 1 (2) 

X» tx 1 - x»/x(0) (X u a lu^t- 1 (3) 

Although Eqs. and (J2J arc applicable for dielectric responses at high frequencies, at 
low frequencies the peak of the response is not often observable; the relaxation process 
is not finalized yet and measured data disregard Eq. © for those processes. Eq. 
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Figure 1. Anomalous low frequency dispersion which is composed of two 
responses, Xi,2i as in Eq. and a dc contribution, a/u>. The thick solid 

line ( ) is the total response. The dashed ( ) and broken ( ) lines 

are the relaxation polarizations. The chain line ( — ■ — ) is the dc conduction 
contribution. 



describes the "normal" low frequency dispersion or an empirical dielectric function, 
also known as constant-phase element (Macdonald 1987). 

X = AO) " 1 and a < 1 (4) 

At the lowest possible frequencies, it is often observed that there is a transition of 
the exponent, a/3 or a, from a low value at high frequencies to a higher value at low 
frequencies beyond a characteristic frequency lu c is passed (Jonscher 1983, Dissado 
and Hill 1984). This phenomenon is called the 'anomalous' low frequency dispersion. 

One can construct an anomalous low frequency dispersion as the summation of 
two dielectric relaxations with t\ < T2 as in the form of Eq. and furthermore, 
a dc conduction contribution a/to can as well be included, where a is normalized 
conductivity-divided by the permittivity of the free space, eo = 1/367T nF/m. 

( s _ Xi(0) , X2(0) Q_ {v s 

Xa[ ' [1 + OwTiH^ [1 + (tWTa) *]* U 

The resulting dielectric response is illustrated in Figure^] Now we discuss the different 
marked-regions in the figure. The real and imaginary parts of the total susceptibilities 
Xa( w ) an d Xai 1 - ) have the same slopes, but different amplitudes on the frequency 
regions above, u> 3> u) c , and below, uj 3> cj c , the critical frequency u) c . In the figure, 
we, therefore, divide the frequency region into three, R±, R2 and R3, which may 
represent a typical experimental data-set. When R3 is considered, the real part of the 
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response, x'ai 1 ^) ~ Xi( w )i is higher than the imaginary part, ~ which 

indicates that 



- l) 2 + I) 2 > 0.5 (6) 

This response in i?3 can be confused with a response containing relaxation polarization 
and dc conduction contributions as expressed below since the finalization of the 
polarization is not clear within this finite frequency range. 

[1 + (ILOT) 01 ^ IU> 

In region i?2, Xa ( w ) is higher than x'ai 1 ^) which shows that in Eq. (JSJ) 



V(« 2 - l) 2 + (fo - I) 2 < 0.5 (8) 

The regions R3 and R2 can separately be described by low frequency disper- 
sions (Jonscher 1983). Consequently, the region R3UR2 is the anomalous low frequency 
dispersion (Dissado and Hill 1984) which is then expressed by summations two power 
laws; 

X (w) = A(iw) 1 - 1 +fl(jw)- i (9) 

where a and b are related to a and j3 as in Eq. l[TJl. and they satisfy the conditions 
in Eqs. JSJ and © in regions i?2 and i?3. A and B depend on the concentration of 
charge carriers, and they are thermally activated constants (Macdonald 1987). The 
intersection of Xa( w ) and Xa( w ) i s the critical frequency, lo c , and it depends on the 
response parameters. 

In region R\ , the presence of the dc conductivity could be observed significantly, 
but, still Xa ( w ) 7^ uuj^ 1 . The dc conductivity contribution hinders the loss peak of 
X2 The bending of or convergence of the real part of the dielectric susceptibility, 
Xa(w)) indicates that the relaxation polarization is to finalize soon. 

In the literature, the low frequency dispersion is associated with materials 
showing defect dominated conduction (hopping conduction) mechanism. Therefore, 
the concentration of defects is critical in the sense that local field changes in the 
materials influences the electrical properties, such as, conductivity of materials that 
do not obey Ohm's law (Sorbello 1984) in the presence of the interaction of charge 
carriers with themselves and with the defects. The micro-structure of the considered 
medium affects the conductivity of the system depending on a length scale which 
is defined as the distance between the charge scatterer or defects (Dissado and 
Hill 1984, Sorbello 1984). Moreover, similar length scales can also be defined as the 
variation of the dielectric permittivity in a continuous random dielectric (Mendelson 
and Schwacher 1984) and the variation of charge density or the electric field. All these 
considerations are valid when there exists some kind of disorder (randomness) in the 
material, i.e., distribution of defects, activation energies, etc. 



3. Modelling considerations 

The electrical properties of the phases, i.e., dielectric permittivity and ohmic 
conductivity are frequency independent, £i,2(k>) = £1,2 and o\_i — constant, in the 
calculations. The phases are linear materials, and the dielectric permittivities, e\ and 
62, of the first (1) and the second (1) phases are taken as 2 and 10, respectively. The 
ohmic conductivities of the media are lOOpS/m ~ 11.3 £0 and IpS/m ~ 1.13 £q. The 
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(a) (b) 

Figure 2. (a) Square and (b) honeycomb 16 X 16 random structures corresponding 
to the same set of lattice computer topology, {%,]). The dark and white polygons 

are phases with different electrical properties. The outer lines ( ) are chosen 

such that the structures are the computation domains. 



considered values of the conductivity and dielectric permittivity yield a polarization 
relaxation around 1 Hz. The concentration of the phases are fifty-fifty qi 2 = 0.5. 

Two-dimensional networks of square and triangular topological arrangements arc 
considered. The networks are composed of 16 x 16 cells, as presented in Figure [3 
and for gi 2 = 0.5, the number of possible different topology configurations is 
- 5.7687 x 10 75 . 

Disordered structures are generated on these lattices by randomly assigned 
material properties such that the distribution of phases are the same for both lattice 
arrangements. This is achieved by generating distributions in the i-j plane, where i 
and j are integers between and 15. The mapping (x-y positions) of the cell centers 
in the square network arc, then, expressed as 

X = z(±+i) y = t(l+j) (10) 



2 

where £ is the length of lattice spacing. In brief, Eq. I|10f> defines the physical 
coordinates (x, y) using the square lattice computer topology, (i,j)- In a similar way 
since the computer topology of the random structures are assumed to be the same, as 
illustrated in Figure the values of is used in the triangular lattice topology to 
assign the physical coordinates (x,y) (Chopard and Droz 1998), 

Evenj: x = £ Q + i\ y = ^fl +j \ (H) 

Oddj: z = £(l + i) y = ^(\+j) ( 12 ) 

(13) 

The inclusions are distributed using a uniform random generator. In this approach, 
each square in the square-network and each honeycomb in the triangular-network are 
treated as the smallest parts of the constituents. Since the corners of squares and 
honeycombs introduce problems in the fem calculations, they are rounded without 
affecting the average shape and area of the inclusions. The rounding process also 
improves the mesh of the computational domain in the finite clement method. We 
use only 1024 random structures which is only ~ 1.7751 x 10~ 73 part of the possible 
topology configurations. 
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Figure 3. Relation between (a) permittivities at high frequencies (uj = 200-7T Hz), 
e, and (b) normalized conductivities at low frequencies (u = 2ir mHz), a for the 
disordered structures with different networks. Subscripts h and s represent the 
honeycomb and square tiles. The open (o) and filled (•) symbols represent the 
considered cases when t\ < 62 and o\ < o"2 and when e\ < £2 and a\ > 02, 
respectively. 



We do not consider the finite-size effects due to the meshing complexity in the 
fem procedure. The first feature of applying two different networks is related to 
the site-bond percolation in the considered networks. When (71.2 = 0.5, the area 
fraction of square network is under the site-percolation threshold (~ 0.592) and 
equal to the bond percolation (0.5). At the same concentration, the area fraction, 
on the other hand, is equal to the site-percolation threshold (0.5) and over the 
bond percolation (~ 0.347) (Stauffer and Aharony 1991). The second feature of our 
geometrical modelling includes the electromagnetic field theory which is not taken 
into consideration in the percolation theory, such as, the conductivity calculations of 
random structures and networks (Stauffer and Aharony 1991). 

4. Results and Discussion 

Dielectric permittivities at constant frequencies 

First, the influence of topological differences are taken into consideration by 
investigating the dielectric permittivities at high (u) = 200-7T Hz) and at low (u) = 
2ir mHz) frequencies. Both analytical mixture expressions and calculations by 
assuming non-percolating regular lattice structures yield steady permittivity values, 
e and e + x(0), at these frequencies (Tuncer et al. 2001a). Two different cases are 
considered in which the permittivities of the phases are kept constant, £\ = 2 and 
£2 = 10. In the first case (case I which is a match composite), the conductivities 
are chosen as o\ = 10~ 12 S/m and er 2 = 10~ 10 S/m. In the other case (CASE II 
which is a reciprocal composite), the conductivities of the phases are interchanged, 
cri = 10- 10 S/m and a 2 = 10~ 12 S/m. 

The results for CASE I and CASE II are illustrated with open (o) and filled 
(•) symbols in Figure |21 In Figure permittivity values obtained from triangular 
lattice arrangement eh are plotted against the permittivities obtained for the square 
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Figure 4. Relation between permittivities, £ = e + x(0) &t lu = 2n mHz for 
the disordered structures with different networks, (a) when ei < €2 and <ri < a"2 
(o) and (b) when ei < t2 and a\ > <T2 (•). Subscripts h and s represent the 
honeycomb and square tiles. 



lattice arrangement e s . There is no linear cross-correlation between the calculated 
permittivities (o symbols in the figure) for the two different topologies for match 
composites (case I). This sort of behavior has also been reported elsewhere for 
simulations performed on 4 x 4 binary mixtures (Tuncer et al. 2001b). However, 
the scattered full symbols (•) for reciprocal composites (CASE II) indicate that there 
is a tendency of correlation since the distribution is like an ellipse. The ellipse is 
slightly parallel to £h = £s diagonal relation but shifted down showing high frequency 
permittivity values for square lattice arrangement has higher polarization compared 
to the triangular arrangement. The regions of > 90 % of the observed values are also 
plotted in the figure as dashed circle and ellipse to display the differences between the 
cases. 

In Figure 03 the normalized conductance values a at 2ir mHz are shown. It 
is important to remind you that at this frequency value non-percolating regular 
structures exhibit clear ohmic losses. For the match composites (case I), the 
distribution of versus a s is narrow compared to CASE II. It is also remarkable 
that there is a clear boundary between match and reciprocal composites such that 
the boundary line is — cr s . This characteristic coincides with that of Tuncer et al. 
(2001b) whose results illustrate that the topological differences can be obtained from 
electrical properties. 

At low frequencies the differences between the dielectric permittivities of the 
considered cases are more significant. In Figures0]i and^t, the dielectric permittivities 
at 2tt mHz are presented for match and reciprocal composites, respectively. In 
Figure^ the influence of structural differences leads higher low frequency permittivity 
values, e — e + x(0), for the square lattice arrangement than triangular lattice with 
honeycomb tiles. When the conductivities of the constituents are interchanged, the 
distribution of versus e s becomes upside-down compared to the CASE I as presented 
in Figure^). The obtained dielectric strengths x(0) are much higher for the reciprocal 
composites 15 < x(0) < 65 than the match composites 1 < x(0) < 9. It is worth to 
mention that the biased distributions in Figure^indicate lower and upper bounds from 
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the structural differences between two lattices; an upper bound for match composites 
and a lower bound for reciprocal composites, respectively. 

4-2. Frequency dependent dielectric properties 

The frequency dependent complex dielectric susceptibilities of two selected structures 
with different lattices are presented in Figure^ The electrical properties of the phases 
are chosen such that e% < €2 and o\ > 02 similar to those reciprocal composites in 
the previous section. We mark regions in which the relaxation changes its behavior as 
in Figure n The Figures |3Jd and|S|i are enlargements of the region R2 in sub-figures 
(a) and (c). The calculated dielectric data do not show the presence of a region 
(i?3 in Figure with low frequency dispersion satisfying the condition in Eq. |j3J) . 
However, we include another region at high frequencies, R4 which shows initiation of 
relaxation processes. At this region both structures have similar dielectric relaxations 
with °c u 2 an d x"( w ) uo ^ 1 indicating the relaxations are Debye-type as 

shown in Figure 01 and|SJ:. However, the dispersions in the intermediate region, R2, 
which is identical to the region presented in Figure ^ can be expressed as distribution 
of relaxation times (Tuncer and Gubahski 2001a, Tuncer et al. 2002a). There is only 
slight differences between two lattice structures at high frequency end of the data 
in i? 2l but at the low frequency side both the structural effects are distinguishable 
from the real x' and the imaginary x" parts of the susceptibilities. The distributions 
presented in Figure^ and0]are once more significant such that at the lowest possible 
frequencies (the region i?i), the structures with honeycomb tiles have higher x'i on 
the other hand, lower x" than square tiles. The structural changes influence the range 
of the low frequency dispersion observed in R2 . The dc conduction contributions are 
significant only in R\ region. 

The dielectric response in Figure |S] has a power law behavior, x(u;) oc (iuj)~ a 
with a = [0.5,0.55] in a short frequency region, 0.1 Hz < u> < 0.3 Hz. When 
a = 0.5 the electrical process is associated with the semi-finite or finite diffusion 
process (Macdonald 1987). Therefore, the mws polarization lead a low frequency 
dispersion. As a result a condition for the 'normal' low frequency dispersion is created 
here with a slightly reciprocal composite. 

To have the "anomalous" low frequency dispersion, the parameters of phase (1) 
are kept constant, ex = 2 and <7\ = 10~ 10 S/m, and the parameters of phase (2) 
are, respectively, chosen (i) 62 = 2000 and 02 = 10 -14 S/m, and (ii) €2 = 20000 
and <72 = 10~ 15 S/m which is an extreme reciprocal composite. The simulated 
dielectric susceptibilities are presented in Figures EJa, and [fjjj. In contrast to the 
previous assumption in Figure [5] the R2 region in this extreme reciprocal composite 
increases in size even more as displayed in Figure This leads to unobservable ohmic 
losses which are present at lower frequencies. Moreover, the flattening of the real part 
(the finalization of polarizations) in Ri is not observable in the considered frequency 
range in the similar region as in Figure^ and[3;. As expected, the x(0) and e values 
are higher than the previous calculations. At the same time, again the power law 
behavior in the region R2 is in xi 10 ) (f-^)~ a form however with a slightly higher 
exponent, a — [0.5,0.6], and satisfying the condition in Equation JHJ). Yet, no sign 
of region R3 is noticeable to the naked eye which would indicate the anomalous low 
frequency dispersion. To show the presence of the anomalous behavior, a parameter, 
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Figure 5. (a) and (b) Dielectric responses obtained for two chosen random 
structures with square (a — •) and honeycomb (□ — ■) tiles. The regions presented 
as i?2 are enlarged in (b) and (d) to emphasize the presence of the low frequency 
dispersions. The open (o — □) and filled (• — I) symbols represent the real, x' 
and imaginary, \" parts of the calculated dielectric susceptibility. 



A, is assigned, 



log 



(14) 



This parameter converges to zero as we approach the critical frequency oj c in Figure^ 
The logarithm of A versus frequency u) is shown in Figures 0: andlSJi for the dielectric 
responses at FiguresEJi andEJi, respectively. The first local minimum at high frequency 
is the first cross-over of \' and x" after the start of the polarization. The region 
between the local minimums in Figures |IJ and 01 indicate frequency interval R3 in 
Figure n m which the condition © is satisfied. The critical frequencies, lj c , for the 
anomalous low frequency dispersions are marked with arrows in the figures. Although, 
the change of the exponent of the power law is not large before and after lu c , as stated 
in the literature (Dissado and Hill 1984), including (i) extra phases in considered 
structures or (ii) using phases with low frequency dispersions or (iii) using a higher 
computation size than 16 x 16 might alter these exponents. 
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Figure 6. Dielectric responses of random structures by keeping phase parameters 
of phase (1) constant, e± = 2 and <j\ = 10~ 10 S/m, and assuming phase (2) 
parameters as (a) t2 = 2000 and ai = 10 — 14 S/m, and (b) 62 = 20000 and 
<T2 = 10 — 15 S/m. The open (o — □) and filled symbols (• — I) represent the real, 
x'(ui) : and imaginary, parts of the dielectric susceptibilities, (c) and (d) 

are plots of the ratio of real, \' (lo), and imaginary, \'' '( w )> P ar ts of the dielectric 
susceptibilities for the considerations in sub-figures (a) and (b), respectively. The 
open circle (o) and open box (□) in (c) and (d) represent the A values obtained 
for square and honeycomb tiles, respectively. 



The influence of topology on the observed low frequency dispersion behavior in 
the middle part of the responses is not that significant, however, on both sides of the 
frequency window, the responses are influenced by the considered lattice structures. 
Finally, since the assumed systems in our simulations are small in size (16 x 16) and 
have the finite smallest tile sizes, a dielectric relaxation of Debye-like behavior with a 
single-time constant is present at high frequencies, the i?4 region. This characteristics 
has also discussed by Macdonald (1987) as the high- and low-frequency extreme limits 
of responses. It has been also presented in Tunccr et al. (2002a), in which the wide 
distribution of relaxation times should start with a single-time constant dielectric 
relaxation (Tuncer et al. 2002a). Above both frequency extremes the low frequency 
dispersion behavior fails. The experimental evidence of this kind of responses in real 
materials, similar to the dielectric relaxations presented in FigureE^, and^s have been 
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Figure 7. Dielectric responses of silicone rubber composites (a) with Aluminum 
Trihydrate powder (b) and glass-beads as inclusions. The open (o — □) and 
filled symbols (• — I) represent the real, x'(uj), and imaginary, x"(u), parts of 
the dielectric susceptibilities. In sub-figure (a) x(ui) of materials with the same 
formulation are presented. In sub-figure (b), the square (□ — I) and circle (o — •) 
symbols represent the the materials with small and large glass-beads, respectively. 



reported in the literature (Jonscher 1983, Dissado and Hill 1984). 
4-3. Dielectric response of filled systems 

In Figure [7] dielectric susceptibilities, x( w )i °f silicone polymer based composite 
materials are presented. In the sub-figures, (a) samples prepared with 50 weight 
% Aluminum Trihydrate powder which has arbitrary-shapes, and (b) with glass-beads 
which are spheroidal are displayed. When the samples with Aluminum Trihydrate 
powder are considered the dielectric responses, presented in Figure [TJi, are similar 
to those presented in Figure [5] of the previous section. In this sample, therefore, the 
conductivity of the filler material (phase 2) might have been considered to have a lower 
conductivity value than the matrix phase (phase 1), 02 < o\, and the permittivity of 
the filler has a higher value than the silicone rubber, t% > ex, a reciprocal composite. 
The data from Tuncer and Gubahski (2000) show that both these conditions are 
true, however, the ratio of permittivities ei/e2 is not as high as we consider in this 
paper. The ratio of the conductivities of the phases is, on the other hand, not an 
easy task to assign. However, the data obtained from different composites with varied 
concentrations of filler indicate that for some concentration value (less than 50 weight 
%) conductivity of the composite decreases. This proves that the conductivity of filler 
<T2 is lower than the matrix when effective medium theories are applied (Sihvola 1999). 
Therefore, the observed low frequency dispersions in Aluminum Trihydrate filled 
silicone rubber samples could be from the MWS polarization. 

In Figure 0d, two responses, in which the composites contain two different 
sizes of glass-beads are illustrated. The larger glass-beads have higher conductivity 
compared with the small glass-beads. We presume that the conductivity of the 
polymer is lower than the conductivity of the large glass-bead, and higher than the 
small glass-beads from the dielectric spectroscopy measurements performed on the 
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polymer. The material with small glass-beads shows an anomalous low frequency 
dispersion (a reciprocal composite), however, the material with large glass-beads 
indicates the presence of two dielectric relaxations and a dc conductivity in the complex 
susceptibility (a match composite). These experimental results together with the 
numerical simulations demonstrate that low frequency dispersion behavior observed 
in dielectric mixtures can be also modeled using computer simulations on disordered 
structures just by assuming the MWS polarization-none intrinsic polarization of 
constituents. 

More illustrative examples in which higher dielectric constants have observed 
are porous materials that are saturated with conductive phases such as salts and 
water (Nettelblad 1996, Papathanassiou and Grammatikakis 2000, Papathanassiou 
2001, Capaccioli et al. 2001, Shen et al. 1990). In these systems real part of the 
permittivity are reported to reach even 10 7 levels as in Nettelblad (1996), however 
in those system the liquid conductivity was responsible for the larger polarization 
observed in the system. In the present work the generated topologies were similar to 
cross-section of porous materials and by assigning material parameters which satisfied 
the condition for "reciprocal mixture", large dielectric permittivity values can be 
obtained as in liquid-solid mixtures. 

5. Conclusions 

The performed fem calculations on binary disordered dielectric mixtures indicate that 
the topology of the system is as important as the intrinsic material properties of the 
constituents to characterize the system. In the simulations 16 x 16 tiles arc generated 
on square and triangular lattices with the same tile arrangement. 

When a match composite, ei < £2 and o\ < o~2, is taken into consideration, 
the topological differences point out no significant relation between the considered 
topological arrangements for permittivity values at high frequencies. This is due to the 
low value of the ratio of the intrinsic permittivity values of the phases. For reciprocal 
composites, ei < 62 and o\ > 0-2, the structural differences does not influence the 
effective dielectric permittivities, e, of the mixture at higher frequencies than the 
characteristic intcrfacial polarization is observed. The relation between e-values are 
nearly-linear. However, other calculated parameters, permittivity at low frequencies 
and dc conductivity, show influence of the structural differences for both match and 
reciprocal composites. This arises from the considered conductivity values which yield 
a high ratio of conductivities. The conductivity also affects the interfacial polarization 
which is visible at very low frequencies. This result confirms that permittivity and 
conductivity can be used as a tool to get structural information confirming the results 
of Tuncer et al. (2001b). At frequencies much lower than the characteristic interfacial 
polarization is observed, the topological differences and the conductivity of the phases 
produce significant differences in the permittivities, e + x(0)- Depending on whether 
the conductivity of the second phase, 172, is higher or lower than ai, the dielectric 
strength of the interfacial polarization, x(0)> differs for the two structures, in which 
in one case, o~\ > 02, the structures with honeycomb tiles have higher polarizations, 
and the opposite is true for structures with square tiles when o~\ < 02 . This indicated 
the importance of the connectivity in the mixture. 

Finally, by altering the conductivity and permittivity of the second phase while 
keeping the electrical properties of the first phase constant, we show that normal 
and anomalous low frequency dispersions can be obtained in dielectric mixtures 
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with frequency independent permittivity and constant dc conductivities. The MWS 
polarization in these systems might be confused with those of normal and anomalous 
low frequency dispersions. The normal low frequency dispersion is observed for 
01 ^> &2 and ei < e 2 - In a like manner, the anomalous low frequency dispersion 
is obtained when o\ 3> 02 and t\ <C £2- Dielectric responses of several polymeric 
composite systems obtained from experiments and of porous materials are compared 
with the responses calculated by the fem. The observed low frequency dispersions 
in those systems might have been due to the MWS polarization and micro-structural 
disorder in the materials which create long paths for charge carriers which are observed 
as to polarization in the dielectric response. 
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